Enhanced transmission versus localization of a light pulse 
by a subwavelength metal slit 



S. V. Kukhlevsky", M. Mechler^ L. Csap6=, K. Janssens'', O. Samek^ 
'^Institute of Physics, University of Pecs, Ifjusdg u. 6, Pecs 7624, Hungary 
South-Trans-Danubian Cooperative Research Centre, 
University of Pecs, Ifjusdg u. 6, Pecs 7624, Hungary 
'^Institute of Mathematics and Information, University of Pecs, Ifjusdg u. 6, Pecs 7624, Hungary 
"^Department of Chemistry, University of Antwerp, Universiteitsplem 1, B-2610 Antwerp, Belgium 
Institute of Spectrochemistry and Applied Spectroscopy, 
Bunsen-Kirchhoff-Str. 11, D-44I39 Dortmund, Germany 

The existence of resonant enhanced transmission and collimation of light waves by subwavelength 
slits in metal films [for example, see T.W. Ebbesen et al.. Nature (London) 391, 667 (1998) and H.J. 
Lezec et al.. Science, 297, 820 (2002)] leads to the basic question: Can a light pulse be enhanced 
and simultaneously localized in space and time by a subwavelength slit? To address this question, 
the spatial distribution of the energy flux of an ultrashort (femtosecond) wave-packet diffracted by 
a subwavelength (nanometer-size) slit was analyzed by using the conventional approach based on 
the Neerhoff and Mur solution of Maxwell's equations. The results show that a light pulse can be 
enhanced by orders of magnitude and simultaneously localized in the near-field difi^raction zone at 
the nm- and fs-scales. Possible applications in nanophotonics are discussed. 



I. INTRODUCTION 



In the last decade nanostructurcd optical elements based on scattering of light waves by subwavelength-size metal ob- 
jects, such as particles and screen holes, have been investigated intensively. The most impressive features of the optical 
elements are resonant enhancement and spatial localization of optical fields by the excitation of electron waves in the 

metal (for example, see the Studieali^i3,4i5,6/,8,9,10,11,12^13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32^33,34,35,36,37)^ 

Recently some nanostructures, namely a single subwavelength slit, a grating with subwavevelength slits, and a sub- 
wavelength slit surrounded by parallel deep and narrow grooves attracted a particular attention of researchers. The 
study of resonant enhanced transmission and collimation of waves in close proximity to a single subwavelength slit 
acting as a microscope probei^ was connected with developing near-field scanning microwave and optical micro- 
scopes with subwavelength resolution^i^*S*i. The resonant transmission of light by a grating with subwavelength slits 
or a subwavelength slit surrounded by grooves is an important effect for nanophotonics^iSiifliii. The transmissivity, 
on the resonance, can be orders of magnitude greater than out of the resonance. It was understood that the enhance- 
ment effect has a two-fold origin: First, the field increases due to a pure geometrical reason, the coupling of incident 
plane waves with waveguide mode resonances located in the slit, and further enhancement arises due to excitation 
of coupled surface plasmon polaritons localized on both surfaces of the slit (grating)i2iii4i^. A dominant mechanism 
responsible for the extraordinary transmission is the resonant excitation of the waveguide mode in the slit giving a 
Fabry-Perot like behaviour—. In addition to the extraordinary transmission, a series of parallel grooves surrounding 
a nanometer-size slit can produce a micrometer-size beam that spreads to an angle of only few degreesSi. The light 
collimation, in this case, is achieved by the excitation of coupled surface plasmon polaritons in the grooves^^. At 
appropriate conditions, a single subwavelength slit flanked by a finite array of grooves can act as a "lens" focusing 
light It should be noted that the diffractive spreading of a beam can be reduced also by using a structured aperture 
or an effective nanolens formed by self-similar linear chain of metal nanospheres^'*'^^. 

New aspects of the problem of resonantly enhanced transmission and collimation of light are revealed when the 
nanostructures are illuminated by an ultra-short (femtosecond) light puls9iSiiLi2ii242&. For instance, in the study—, 
the unique possibility of concentrating the energy of an ultrafast excitation of an "engineered" or a random nanosystem 
in a small part of the whole system by means of phase modulation of the exciting fs-pulse was predicted. The study- 
theoretically demonstrated the feasibility of nm-scale localization and distortion-free transmission of fs visible pulses 
by a single metal slit, and further suggested the feasibility of simultaneous super resolution in space and time of the 
near-field scanning optical microscopy (NSOM). The quasi-diffraction-free optics based on transmission of pulses by a 
subwavelength nano-slit has been suggested to extend the operation principle of a 2D NSOM to the "not-too-distant" 
field regime (up to ~ 0.5 wavelength)^^. Some interesting effects, namely the pulse delay and long living resonant 
excitations of electromagnetic fields in the resonant-transition gratings were recently described in the studiesi^iSS. 

The great interest to resonant enhanced transmission, spatial localization (collimation) of continuous waves and 
light pulses by subwavelength metal slits leads to the basic question: Can a light pulse be enhanced and simultaneously 
localized in space and time by a subwavelength slit? If the field enhancement can be achieved together with nm-scale 
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spatial and fs-scale temporal localizations, this could greatly increase the potentials of the nanoslit systems in high- 
resolution applications, especially in near-field scanning microscopy and spectroscopy. In the present article we test 
whether the resonant enhancement could only be obtained at the expense of the spatial and temporal broadening of 
a light wavepacket. To address this question, the spatial distribution of the energy flux of an ultrashort (fs) pulse 
diffracted by a subwavelength (nanosizcd) slit in a thick metal film of perfect conductivity will be analyzed by using 
the conventional approach based on the Neerhoff and Mur solution of Maxwell's equations. In short, we first will 
summarize the theoretical development of Neerhoff and Mur (Section II) and the model will then be used to calculate 
the spatial distribution of the energy fiux of the transmitted pulse (wavepacket) under various regimes of the near-field 
diffraction (Section III). We will show that a light pulse can be enhanced by orders of magnitude and simultaneously 
localized in the near-field diffraction zone at the nm- and fs-scales. The implications of the results for diffraction- 
unlimited near- and far-field optics will then be discussed. In Section IV we summarize the results and present the 
conclusions. 

II. THEORETICAL BACKGROUND 

An adequate description of the transmission of light through a subwavelength nano-sizcd slit in a thick metal film 
requires solution of Maxwell's equations with complicated boundary conditions. The light-slit interaction problem 
even for a continuous wave can be solved only by extended two-dimensional {x, z) numerical computations. The 
three-dimensional {x,z,t) character of the pulse-slit interaction makes the numerical analysis even more difficult. 
Let us consider the near-field diffraction of an ultrashort pulse (wave-packet) by a subwavelength slit in a thick 
metal screen of perfect conductivity by using the conventional approach based on the Neerhoff and Mur solution 
of Maxwell's equations. Before presenting a treatment of the problem for a wave-packet, we briefly describe the 
Neerhoff and Mur formulation-''*3i for a continuous wave (a Fourier w-component of a wavepacket). The transmission 
of a plane continuous wave through a slit (waveguide) of width 2a in a screen of thickness b is considered. The 
perfectly conducting nonmagnetic screen placed in vacuum is illuminated by a normally incident plane wave under 
TM polarization (magnetic-field vector parallel to the slit), as shown in Fig. 2] The magnetic field of the wave is 
assumed to be time harmonic and constant in the y direction: 

H{x, y, z, t) — U{x, z)ejip{—iujt)ey. (1) 
The electric field of the wave is found from the scalar field U{x, z) using Maxwell's equations: 

Ex{x, z, t) — dzU{x, z)exp(— iwt), (2) 

Ey{x,z,t) = 0. (3) 
tc 

Ez{x,z,t) — dxU{x, z)exp{—iujt). (4) 

Notice that the restrictions in Eq. ^reduce the diffraction problem to one involving a single scalar field U{x, z) in two 
dimensions. The field is represented by C/j(x, z) (j'=l,2,3 in each of the three regions I, II and III). The field satisfies 
the Helmholtz equation: 

(V^ + k])Uj = 0, (5) 
where j = 1, 2, 3. In region I, the field Ui{x, z) is decomposed into three components: 

Ui{x, z) = U\x, z) + Wix, z) ^ V^ix, z), (6) 

each of which satisfies the Helmholtz equation, f/* represents the incident field, which is assumed to be a plane wave 
of unit amplitude: 

U^{x,z) — cx]i{-~ikiz). (7) 
[/'' denotes the field that would be reflected if there were no slit in the screen and thus satisfies 



Vix^z) = W{x,2b-z). 



(8) 
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describes the diffracted field in region I due to the presence of the slit. With the above set of equations and 
standard boundary conditions for a perfectly conducting screen, a unique solution exists for the diffraction problem. 
To find the field, the 2-dimensional Green's theorem is applied with one function given by U{x^ z) and the other by a 
conventional Green's function: 

{\I^ + k])Gj^^5{x-x',z-z'), (9) 

where (x, z) refers to a field point of interest; x', z' are integration variables, j = 1,2,3. Since Uj satisfies the Helmholtz 
equation, Green's theorem reduces to 

U{x,z)^ j [GdnU -UdnG)dS. (10) 

J Boundary 

The unknown fields J7''(x, z), U^{x^z) and U2{x,z) are found using the reduced Green's theorem and boundary 
conditions on G 

U'^{x,z)^-—[ Gi{x,z-x',b)DUb{x')dx' (11) 

£2 J-a 



for b < z < oo, 



for —oo < z < 0, 



Uiix.z)^— I G3{x,z;x',0)DUo{x')dx' (12) 

£2 



U2{x, z) = - [G2(x, z- x\ 0)DUo{x') ~ Uo{x')d,,G2{x, z; x\ z')\,^„+]dx' 

J —a 

+ f [G2ix,z;x',b)DUbix')-Ub{x')dz'G2{x,z;x\z%^t-]dx' (13) 



for < a and < z < b. The boundary fields in Eqs. II 11131 are defined by 



Uo{x) = U2{x,z)U^o+, (14) 

DUoix) ^d,U2ix,z)\,^o+, (15) 

Ui,{x)^U2{x,z)\,^,-, (16) 

DUb{x)^d,U2{x,z)\,^h-. (17) 



In regions I and III the two Green's functions in Eqs. II II and [T^ are given by 



Gi{x,z;x',z') = '-[H^^''\hR) + H^,'\k,R% (18) 
G3(x,z;x',z') = '-[4'\hR) + H^'\hR")], (19) 



with 



R=[{x-x')^ + {z- z')^]^^^, (20) 
R' ^ [{x - x'f + (z + z' - 2bf]^^\ (21) 
R" ^[{x-x'f + iz + z'f]'/^, (22) 



where Hq^^ is the Hankel function. In region II, the Green's function in Eq. ^|is given by 



G2{x,z;x',z') = exp(i7o|z- z'l) + -^Y] 7™^ 
4a7o 20^^™=! 

xcos[77i7r(a: + a)/2a]cos[m7r(a;' + a) /2a] 

xexp(i7,„|0 - z'l), (23) 
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where 7™ = [A;| - (m7r/2a)2]i/2_ The field can be found at any point once the four unknown functions in Eas. 1141 Ffl 
have been determined. The functions are completely determined by a set of four integral equations: 



2Ul{x) - Ub{x) = — / Gi{x,b;x',b)DUt{x')dx' 

f-2 



Uoix) 



£3 

£2 



G3{x,0;x',0)DUoix')dx' 



(24) 
(25) 



[G2{x, b; x', 0)DUo{x') - Uo{x')d,>G2{x, 6; x' , z')\,^o^]dx' 

[G2(x, b; x , b)DUb{x')]dx' , 



(26) 



-Uoix) 



[G2{x, 0; x', b)DUb{x') - Ub{x')d,,G2{x, 0; x', z')U^b-]da;' 

[G2{x,0]x' ,0)DUo{x')]dx' , 



where \x\ < a, and 



Ul{x) = exp{~ikib). 



(27) 



(28) 



The coupled integral equations I24I27I for the four boundary functions are solved numerically. The magnetic H{x, z, t) 
and electric E{x,z,t) fields of the diffracted wave in region III are found by using Eq. El The fields are given by 



#(x,Z,i) = *T7-E^0 



N 



(1) 



fcsy {x - XjY + 
X (DUq) jC'K.p{--iujt)ey 



(29) 



E^{x,z,t) 



a V<^3 



N 

E 



N £2 ^ ^(a; - XjY + z2 



(1) 



i?j;(a;, z, i) = 0, 



E^{x,z,t) 



N 

E 



N £2 ^^(a;-x,)2 + z2 



=i7 



(1) 



kz\J {x — XjY + z' 
X (D[/o)jexp(— io^i) 

k^yj {x — Xjf- + z2 
X (DC/o)jexp(— iwt), 



(30) 
(31) 



(32) 



where a;^ = 2a(j — l/2)/N ~ a, j — 1,2, ...,-/V; > 2a/z; iJ;^^-' is the Hankel function. The coefhcients (DUo)j are 
found by solving numerically the four integral equations I24I27I For more details of the model and the numerical 
solution of the coupled integral equations I24I27I see refs>ii2i. 

Let us now consider the diffraction of an ultra-short pulse (wave packet). The magnetic field of the incident pulse 
is assumed to be Gaussian-shaped in time and both polarized and constant in the y direction: 



H{x,y,z,t) — U{x, z)exp[— 2 \n{2){t / T)'^]exp{~iujot)ey , 



(33) 



where t is the pulse duration and luq = 2ttc/ Xq is the central frequency. The pulse can be composed in the wave-packet 
form of a Fourier time expansion (for example, see ref^ii^): 



H{x,y,z,t)= / H{x, z,uj)exp{—iLut)duj. 



(34) 
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The electric and magnetic fields of the diffracted pulse are found by using the expressions H1I32|I for each Fourier 
w-component of the wave-packet (I34II . This algorithm is implemented numerically by using the discrete Fast Fourier- 
Transform (FFT) instead of the integral (j^H). The spectral interval [cOmin,^max] and the sampling points cui are 
optimized by matching the FFT result to the original function (|33|l . 

The above presented approach deals with the incident waves having TM polarization. This polarization is considered 
for the following reasons. According to the theory of waveguides, the vectorial wave equations for this polarization 
can be reduced to one scalar equation describing the magnetic field H of TM modes. The electric component E of 
these modes is found using the field H and Maxwell's equations. The reduction simplifies the diffraction problem to 
one involving a single scalar field in only two dimensions. The TM scalar equation for the component H is decoupled 
from the similar scalar equation describing the field E of TE (transverse electric) modes. Hence, the formalism 
works analogously for TE polarization exchanging the E and H fields. Moreover, in the case of perfectly conducting 
nonmagnetic screen placed in vacuum considered in the present paper, the symmetry of wave equations indicates that 
large transmission coefficients (enhancement effect) do exist at the same experimental conditions for the TM and TE 
polarizations. 



III. NUMERICAL ANALYSIS AND DISCUSSION 



In this section, we test whether a light pulse can be resonantly enhanced and simultaneously localized in space and 
time by a subwavelength nano-sized metal slit. To address this question, the spatial distribution of the energy flux of 
the transmitted pulse under various regimes of the near-field diffraction is analyzed numerically. The electric E and 
magnetic H fields of the transmitted pulse in the near-field diffraction zone are computed by solving the equations 1^ 
I32|l for each Fourier w-component of the wave-packet (j34|l . The amplitude of a FFT w-component of the wave-packet 
transmitted through the slit depends on the wavelength A = l-ncjio. Owing to the dispersion, the Fourier spectra 
of the transmitted wave-packet changes leading to modification of the pulse width and duration. The dispersion 
of a continuous wave is usually described by the normalized transmission coefficient Tc^(A), which is calculated by 
integrating the normalized energy flux Szj S\ over the slit widtbl^:: 



4a 



lim [{E^m 



(35) 



where SI is the energy flux of the incident wave of unit amplitude; Sz is the transmitted flux. Equation 1351 which 
includes the evanescent modes that decay in the far zone, deflncs the transmission coefficient in the near-field zone. 
In order to establish guidelines for the results of our numerical analysis, we computed the transmission coefficient 
Tcu,(A, a, 6) for a continuous wave (Fourier w-component) as a function of screen thickness b and/or wavelength A 
for different values of slit width 2a. Throughout the computations, the magnitude of the incident magnetic field is 
assumed to be equal to 1. We consider a perfectly conducting nonmagnetic screen placed in vacuum (ei = £2 = £3 = 1) 
As an example, the dependence Tew = Tcw{b) computed for the wavelength A = 800 nm and the slit width 2a = 25 nm 
is shown in Fig. |21 The dispersion Tew — TcwW for 2a — 25 nm and different values of the screen thickness b is 
presented in Fig.|3| In Fig.|21 we note the transmission resonances of A/2 periodicity with the peak heights Tc^w A/27ra 
at the resonances. Notice that the resonance positions and the peak heights are in agreement with the result^"^. The 
resonance peaks appear with spacing of 400 nm, but at values somewhat smaller than 400, 800, 1200 nm, etc. In 
order to better understand the shift of the resonant wavelengths from the naively expected values we have derived a 
simple analytical formula H3fci|) for the transmission coefficient T^w = Tcw{ci,b, X) of a narrow slit (the fundamental 
mode only is retained): 



Tcwia, b,X)^^ [(Re {D (6, A)))' + (Im {D {b, X))f 



(36) 



where 



^(^A) = — 



r/A 



A{X) 



27r/A 



(37) 



and 



A{X) la [ff (A) + I (Ho (A) • ^ (A) + Hi (A) • H^^^] (A) 



(38) 
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In Eq. OHl i?o^^(A) and h[^\x) are the Hankel-functions, while Ho(A) and Hi(A) are the Struve-functions. It is inter- 
esting to compare the "near-field" {z << X) transmittance with the conventional "far-zone" (z >> A) transmittance 
measured in experiments. The transmittance in the far-field zone can be described by the following formula: 

/e~ f°° 

Tc^u - / Hm [{E,H* + E*Hy)]dx, (39) 

Figure 121 compares the transmission coefficients Tew — Tcw{b) calculated by using the formula H36|l with the values 
given by the numerical solution of the equations H24l27llS3|) and the values obtained from the evaluation of (|39|l . We 
notice that the results are practically undistinguishable. Analysis of the denominator of the formula (|36|1 indicates 
that for 2a/A small enough, the transmission Tew — Tcw(b) will exhibit the maximums around wavelengths ~ A/2. 
The shifts of the resonance wavelengths from the value A/2 are caused by the wavelength dependent terms in the 
denominator of Eq. 1361 

The dispersion Tew = Tcw{X) presented in Fig. |21 indicates the wave-slit interaction behaviour, which is similar to 
those of a Fabry-Perot resonator. The transmission resonance peaks, however, have a systematic shift towards longer 
wavelengths. Analysis of the denominator of the formula (|36|l indicates that for 2a/ A small enough, the transmission 
Tew = Tew{X) will exhibit the Fabry-Perot like maximums around wavelengths where sin(fc6) is zero. Such a behaviour 
has already been described in refs 4^*2^21. The shifts of the resonance wavelengths from the Fabry-Perot resonances 
are caused by the wavelength dependent terms in the denominator of Eq. 1211 It should be noted that the values 
of the shifts calculated using Eq. 1211 are in good agreement with the shifts calculated using the analytical formula 
(8) of the study — . Our computations showed that the peak heights at the main (strongest) resonant wavelength 
A^ (in the case of Fig. |21 A^= 500 or 800 nm) are given by TcwiXg , a)^X§ /2na. Notice that the Fabry-Perot like 
behaviour of the transmission coefficient is in agreement with analytical and experimental results published earlieriSL^. 
The enhancement coefficient {Tew ~ 30) calculated using Eas. 1351 and 1361 however, is different from the attenuation 
(Tew < 1) predicted by the study ^''^ . The one order difference between the experimental value and the enhancement 
calculated in the present article is attributed to the small transverse dimension D (D ^ X/2) of the metal plates that 
form the slit. In our calculations we studied a slit in infinite transverse dimension {D = oo). 

Analysis of Fig. |2| indicates that a minimum thickness of the screen is required to get the waveguide resonance 
inside the slit at a given wavelength. The result is in agreement with the study which demonstrated that there 
is no transmission peaks at the condition b < A/2. Notice that the transmission enhancement mediated by surface 
plasmons does exist at considerably smaller thicknesses of the metal in comparison with the thickness required for 
the waveguide resonance. The surface plasmons/polaritons are excited in the skin, the depth of which is about the 
extinction length of the energy decay of electromagnetic wave in the metal. For instance, the extinction length in 
an aluminum screen is ~ 4 nm for A = 800 nm. Therefore, for not too narrow slits (2a = 25 nm) in thick screens 
{b = 200 and 350 nm) considered in the paper, the finite skin depth does not fundamentally modify the process of the 
extraordinary optical transmission. At the wavelength 800 nm, the 10-times enhancement (Tcu,(a, A) « A/27ra « 10) 
is limited by the slit width (2a = 25 nm). In the far-infrared region (A ~ 100 ^m) several orders of magnitude 
enhancement can be achieved at the same experimental conditions. It should be noted that the optimal choice of 
parameters has been discussed in the literature and the obtained enhancement by the factors 10 and 1000 are typical 
for continuous waves in the optical and far-infrared spectral ranges. In the case of narrower and thinner slits, however, 
the influence of the finite conductivity effects on the transmission and localization of a pulse should be taken into 
account -^sSi. 

The existence of transmission resonances for Fourier w-components of a wave-packet leads to the question: What 
effect the resonant enhancement has on the spatial and temporal localization of a light pulse? Presumably, the high 
transmission at resonance occurs when the system efficiently channels Fourier-components of the wave-packet from a 
wide area through the slit. At resonance, one might assume that if the energy flow is symmetric about the screen, 
the pulse width and duration should increase very rapidly past the screen. Thus the large pulse strength associated 
with resonance could only be obtained at the expense of the spatial and temporal broadening of the wave-packet. 
To test this hypothesis, the spatial distributions of the energy flux of a transmitted wave-packet were computed for 
different slit thicknesses corresponding to the resonance and anti-resonance position. The characteristic difference 
between the resonant and non-resonant transmissions could be understood better through a single figure of the field 
distributions in all regions I, II and III. However, it seems to be impossible to do this, because the visualization of 
a pulse by a single figure requires a 4-dimensional (U, x, z, t) coordinate system or a great number of 3-dimensional 
figures at different locations and suitable fixed times in all the regions. Therefore we present the eight 3-dimensional 
distributions l^a)-(d) and|Ha)-(d)) only in the zone of main interest (region III). As an example. Figs. 0] (a) and 
121 (a) show the energy flux of the anti-resonantly transmitted pulses. Figures 21(b) and [21(b) correspond to the case 
of the waveguide- mode resonance in the slit. Figures ^ and (21 show the transmitted pulses at the distances \z\ — a/2 
and a, respectively. The rigorous analysis^ showed that the number of modes required for the accurate computation 
of the transmittance, near-field distribution, and other optical characteristics of the system is given by TV > 2a/ z, 
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where N is the number of the waveguide modes and z is the distance from the screen. Therefore at the distances 
\z\ — a/2 and a the calculations required at least 4 and 2 modes, respectively. 

The shape and intensity of an output pulse depends on the slit parameters and the spectral width of the pulse. 
For narrow slits, the spectral width of a 100-fs input pulse is smaller compared to the spectral width of the resonant 
transmission (see, Fig. O. The comparison of the flux distribution presented in Fig. 0] (a) with that of Fig. ^ (b) 
shows that, for the parameter values adopted, a transmitted wavepacket is enhanced by one order of magnitude and 
simultaneously localized in the 25-nm and 100-fs domains of the near-field diffraction zone. The shapes of the intensity 
distributions of the output pulses are very much the same off and on resonance. The figures differ only in the order 
of magnitude of Sz- Thus at the distance \z\ = a/2, the slit resonantly enhances the intensity of the pulse without its 
spatial and temporal broadening. The result can be easily understood by considering the dispersion properties of the 
slit. For the screen thickness h = 200 nm, the amplitudes of the Fourier-components of the wave-packet, whose central 
wavelength Aq is detuned from the main (at 500 nm) resonance, are practically unchanged in the wavelength region 
near 800 nm (see, curve B in Fig.lSJ. This provides the dispersion- and distortion-free non-resonant transmission of 
the wave-packet (Figs.QJa) andE^a)). In the case of the thicker screen (b = 350 nm), the slit transmission experiences 
strong mode-coupling regime at the wavelengths near to 800 nm (see curve A of Fig. Ojl that leads to a profound and 
uniform enhancement of amplitudes of all of the Fourier tj-components of the wave-packet (see curve C in Fig. |2l . 
Thus, the slit resonantly enhances by one order of magnitude the intensity of the pulse without its spatial and temporal 
broadenings (Figs.^fb) and^fb)). Also, notice that at the distance \z\ — a (Figs. EJa) and[SJb)), both the resonantly 
and anti-resonantly transmitted pulses experience natural spatial broadening in the transverse direction, while their 
durations are practically unchanged. The spectral width of a 5-fs input pulse, however, is bigger than the spectral 
width of the resonant transmission (Fig. |3Jl. In this case the durations of the resonantly transmitted 5-fs pulses are 
longer in comparison with the pulses transmitted off the resonance (Figs. EIc), (d) and respectively [S{c), (d)). We 
also notice that the enhancement of the intensity of 5-fs pulses is approximately two times lower in comparison with 
the 100-fs pulses. The temporal broadening of the pulse and the decrease of the enhancement indicates a natural 
limitation for the simultaneous temporal and spatial localization of a pulse together with its enhancement. 

By comparing the data for anti-resonant and resonant transmissions presented in Figs. ^ and [S] one can see that at 
the appropriate values of the distance \z\ and the wave-packet spectral width Aw the resonance effect docs not influence 
the spatial and temporal localization of the wave-packet. To verify this somewhat unexpected result, the FWHMs of 
the transmitted pulse in the transverse and longitudinal directions were calculated for different values of the slit width 
2a, central wavelength Aq and pulse duration rwl/Aojp as a function of screen thickness b at two particular near-field 
distances \z\ = a/2 and a from the screen. It was seen that, at the dispersion- free resonant transmission condition 
Awp < Awr, the transmitted pulse indeed does not experience temporal broadening. Thus the temporal localization 
associated with the duration r of the incident pulse remains practically unchanged under the transmission. The value 
of T is determined by the dispersion- free condition r^l/Awp > l/Aw^, where Atj^ = Awr(a) practically does not 
depend on the screen thickness h. We found that the energy flux of the transmitted wave-packet can be enhanced 
by a factor Tc^(A^, a)«A^/27ra by the appropriate adjusting of the screen thickness 6 — b{X§), for an example see 
Figs. 1351 Thus the wave-packet can be enhanced by a factor A^/27ra and simultaneously localized in the time domain 
at the r — r(a) scale. It was also seen that the FWHM of the transmitted pulse in the transverse direction depends on 
the wave-packet central wavelength Aq and the distance z from the slit. Nevertheless, the FWHM of the transmitted 
pulse can be always reduced to the value 2a by the appropriate decreasing of the distance 1^1 = \z{a)\ from the screen 
(|z| — a/2, in the case of Fig. Thus high transmission can be achieved without concurrent loss in the degree of 
temporal and spatial localizations of the pulse. In retrospect, this result is reasonable, since the symmetry of the 
problem for a time-harmonic continuous wave (Fourier w-component of a wave-packet) is disrupted by the presence 
of the initial and reflected fields in addition to the diffracted field on one side (Eq. |HJ). As the thickness changes, the 
fields and U3 change only in magnitude, but the field Ui changes in distribution as well since it involves the sum 
of with unchanging fields W and U^. At resonance, the distribution of Ui leads to channeling of the radiation, 
but the distribution of t/3 remains unaffected. By the appropriate adjusting of the slit-pulse parameters a light pulse 
can be enhanced by orders of magnitude and simultaneously localized in the near-field diffraction zone at the nm- 
and fs-scales. 

The limitations of the above analysis must be considered before the results are used for a particular experimental 
device. The resonant enhancement with simultaneous nm-scale spatial and fs-scale temporal localizations of a light 
by a subwavelength metal nano-slit is a consequence of the assumption of the screen perfect conductivity. The slit 
can be made of perfectly conductive (at low temperatures) materials. In the context of current technology, however, 
the use of conventional materials like metal films at a room temperature is more practical. Notice that a metal 
can be considered as a perfect conductor in the microwave range. As a general criterion, the perfect conductivity 
assumption should remain valid as far as the slit width and the screen thickness exceed the extinction length for 
the Fourier tj-components of a wavepacket within the metal. The light intensity decays in the metal screen at the 
rate of Ig = Iq exp(— 6/5), where S — S{X) is the extinction length in the screen. The aluminum has the largest 
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opacity {S < 11 nm) in the spectral region A > 100 nm [20]. The extinction length increases from 11 to 220 nm 
with decreasing the wavelength from 100 to 50 nm. Hence, the perfect conductivity is a very good approximation in 
a situation involving a relatively thick (6 > 25 nm) aluminum screen and a wave-packet of the duration rsal/Aojp 
having the Fourier components in the spectral region A > 100 nm. However, in the case of thinner screens, shorter 
pulses and smaller central wavelengths of wave-packets, the metal films are not completely opaque. This would reduce 
a value of spatial localization of a pulse due to passage of the light through the screen in the region away from the 
slit. Moreover, the phase shifts of the Fourier components along the propagation path caused by the skin effect can 
modify the enhancement coefficient and temporal localization properties of the slit. 

The above analysis is directly applicable to the two-dimensional near-field scanning optical microscopy and spec- 
troscopy. Although the computations were performed in the case of normal incidence, the preliminary analysis shows 
that in the case of oblique incidence the enhancement effect can be kept without further spatial and temporal broad- 
ening of the pulse. In a conventional 2D NSOM, a subwavelength (2a < A) slit illuminated by a continuous wave 
is used as a near- field {\z\ < a) light source providing the nm-scale resolution in space^iii^. The non-resonant 
transmission of fs pulses could provide ultrahigh resolution of 2D NSOM simultaneously in space and tim o^Ti^? . The 
above-described resonantly enhanced transmission together with nm- and fs-scale localizations in the space and time 
of a pulse could greatly increase the potentials of the 2D near-field scanning optical microscopy and spectroscopy, 
especially in high-resolution applications. However, as a model for a NSOM-tip a hole (quadratic, rectangular or 
circular) would be more appropriate than a slit-type waveguide. The consequences of the limitation of region II also 
in y-direction will result in faster attenuation of the amplitude of the waveguide modes compared to the slit-type 
waveguide modes—. This should affect the enhancement and spatial and temporal broadening of the pulse. Moreover, 
the walls of the 2D slit have to be parallel over all the thickness to provide the resonance enhancement effect. Pulled 
NSOM 3D-tips, however, usually have conical ends. In the case of the tapered NSOM 3D waveguides one should 
also take into account the decrease of enhancement effect with increasing the waveguide taper. It should be also 
noted that the high transmission (Tcu)(A^, a)R:!A^/27ra) of a pulse can be achieved without concurrent loss in the 
temporal and spatial localizations of the pulse only at short (|z| — \z{a)\) distances from the slit. The presence of 
a microscopic sample (a molecule, for example) placed at the short distance in strong interaction with NSOM slit 
modifies the boundary conditions. In the case of strong slit-sample-pulse interaction, which takes place at the distance 
\z\ << 0.1a, the response function accounting for the modification of the quantum mechanical behaviour of the sample 
should be taken into consideration. The potential applications of the effect of the resonantly enhanced transmission 
together with nm- and fs-scale localizations of a pulse are not limited to near-field microscopy and spectroscopy. 
Broadly speaking, the effect concerns all physical phenomena and photonic applications involving a transmission of 
light by a single subwavelength nano-slit, a grating with subwavevelength slits and a subwavelength slit surrounded 

by parallel grooves (see the Studieali^i3,4,5,6,7,8,9,10,ll,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37 

references therein). For instance, the effect could be used for sensors, communications, optical switching devices and 
microscopes. 

IV. CONCLUSION 

In the present article we have considered the question whether a light pulse can be enhanced and simultaneously 
localized in space and time by a subwavelength metal nano-slit. To address this question, the spatial distributions 
of the energy flux of an ultrashort (fs) pulse diffracted by a subwavelength (nanosized) slit in a thick metal screen of 
perfect conductivity have been analyzed by using the conventional approach based on the Nccrhoff and Mur solution 
of Maxwell's equations. The analysis of the spatial distributions for various regimes of the near-field diffraction 
demonstrated that the energy flux of a wavepacket can be enhanced by orders of magnitude and simultaneously 
localized in the near-field diffraction zone at the nm- and fs-scales. The extraordinary transmission, together with 
nm- and fs-scale localizations of the light pulse, make the nano-slit structures attractive for many photonic purposes, 
such as sensors, communication, optical switching devices and NSOM. We also believe that the addressing of the 
above-mentioned basic question gains insight into the physics of near-field resonant diffraction. 

Acknowledgments 

This study was supported by the Fifth Framework of the European Commission (Financial support from the EC 
for shared-cost RTD actions: research and technological development projects, demonstration projects and combined 
projects. Contract NG6RD-CT-2001-00602). The authors thank the Computing Services Centre, Faculty of Science, 



9 



University of Pecs, for providing computational resources. 



F. L. NeerhofT, G. Mur G, Appl. Sci. Res. 28, 73 (1973). 

R.F. Harrington, D.T. AucklandJ, IEEE Trans. Antennas Propag AP28, 616 (1980). 
E. Betzig, A. Harootunian, A. Lewis, M. Isaacson, Appl. Opt. 25, 1890 (1986). 
E.A. Ash, G. NichoUs, Nature (London) 237, 510 (1972). 
A. Lewis, K. Lieberman, Nature (London) 354, 214 (1991). 

E. Betzig, J.K. Trautman, T.D. Harris, J.S. Weiner, R.L. Kostelak, Science 251, 1468 (1991). 

D. W. Pohl, D. Courjon, Eds., Near Field Optics, NATO ASI Series E Vol. 242 (The Netherlands, Dordrecht: Kluwer, 1993). 
T.W. Ebbesen, H.J. Lezec, H.F. Ghaemi, T. Thio, PA. Wolff, Nature (London) 391, 667 (1998). 

^ H.J. Lezec, A. Degiron, E. Devaux, R.A. Linke, Martin-Moreno, F.J. Garcia- Vidal, T.W. Ebbesen, Science 297, 820 (2002). 

J. A. Porto, F.J. Garcia-Vidal, J.B. Pendry, Phys. Rev. Lett. 83, 2845 (1999). 
" A.P. Hibbings, J.R. Sambles, C.R. Lawrence, Appl. Phys. Lett. 81, 4661 (2002). 

F. J. Garcia-Vidal, H.J. Lezec, T.W. Ebbesen, L. Martin-Moreno, Phys. Rev. Lett. 90, 231901 (2003). 
F.J. Garcia-Vidal, L. Martin- Moreno, H.J. Lezec, T.W. Ebbesen, Appl. Phys. Lett. 83, 4500 (2003). 
A. Dogariu, M. Hsu, L.J. Wang, Opt. Comm. 220, 223 (2003). 

K.R. Li, M.I. Stockman, D.J. Bergman, Phys. Rev. Lett. 91, 227402 (2003). 
" M.I. Stockman, S.V. Faleev, D.J. Bergman, Phys. Rev. Lett. 88, 067402 (2002). 
S.V. Kukhlevsky, M. Mechler, L. Csapo, K. Janssens, Phys. Lett. A 319, 439 (2003). 
S.V. Kukhlevsky, M. Mechler, Opt. Comm. 231, 35 (2004). 
P. N. Stavrinou, Phys. Rev. E 68, 06604 (2003). 
A.G. Borisov, S.V. Shabanov, arXiv:physics/0312103 (2003). 
Y. Takakura, Phys. Rev. Lett. 86, 5601 (2001). 
F.Z. Yang, J.R. Sambles, Phys. Rev. Lett. 89, 063901 (2002). 
U. Schroter, D. Heitmann, Phys. Rev. B 58, 15419 (1998). 
M.M.J. Treacy, Phys. Rev. Lett. 75, 606 (1999). 

L. Salomon, F.D. Grillot, A.V. Zayats, F. de Fornel, Phys. Rev. Lett. 86, 1110 (2001). 
Q. Cao, P. Lalanne, Phys. Rev. Lett. 88, 057403 (2002). 

I.I. Smolyaninov, A.V. Zayats, A. Stanishevsky, C.C. Davis, Phys. Rev. B 66, 205414 (2002). 
A. Barbara, P. Quemerais, E. Bustarret, T. Lopez-Rios, Phys. Rev. B 66, 161403 (2002). 

E. Altewischer, M.P. van Exter, J.P. Woerdman, Nature (London) 418, 304 (2002). 
A.M. Dykhne, A.K. Sarychev, V.M. Shalaev, Phys. Rev. B 67, 195402 (2003). 
J.M. Steele, C.E. Moran, A. Lee, CM. Aguirre, N.J. Halas, Phys. Rev. B 68, 205103 (2003). 
A. Naweed, F. Baumann, W.A. Bailey, A.S. Karakashian, W.D. Goodhue, J. Opt. Soc. Am. B 20, 2534 (2003). 
Z.J. Sun, Y.S. Jung, H.K. Kim, Appl. Phys. Lett. 83, 3021 (2003). 
X.L. Shi, L. Hesselink, R.L. Thornton, Opt. Lett. 28, 1320 (2003). 

H.F. Schouten, T.D. Visser, D. Lenstra, H. Blok, Phys. Rev. E 67, 036608 (2003). 
'■^^ A. Nahata, R.A. Linke, T. Ishi, K. Ohashi, Opt. Lett. 28, 423 (2003). 
3^ G. Gomez-Santos, Phys. Rev. Lett. 90, 077401 (2003). 



10 



14 



10 



H 



Z 

M - 



Incident wave 

H k 
^ *- 



a 



X 



IE 



III 



Diffracted 
wave 



FIG. 1: Propagation of a continuous wave through a subwavelength nano-sized slit in a thick metal film. 
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FIG. 2: The transmission coefficient Tew for a continuous wave (tj-Fourier component of a wave-packet) as a function of screen 
thickness b computed for the wavelength A=800 nm and the slit width 2a — 25 nm. Curve A, B and C correspond to the 
numerical near-field formula H,'i5|l . the analytical near-field formula 1)36^ and the analytical far-field formula Ijiffl^ . respectively. 
For computational reasons, in the case of curve C the limits of the integral were chosen to -lOOOA and lOOOA instead of — oo 
and oo; the computation was performed ai z — — lOA. 
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FIG. 3; The dispersion T^w = Tciii(A) for a continuous wave (oj-Fourier component of a wave-packet) computed for the slit 
width 2a = 25 nm and different values of the screen thickness fo: A - 350 nm and B - 200 nm. The Fourier spectra (curves 
C and D) are presented for the comparison. Curves C and D show the Fourier spectra of incident wave-packets with central 
wavelength Ao — 800 nm and duration r = 100 fs and r = 5 fs, respectively, which were used in the computations presented in 
Fig. Hand H 
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FIG. 4: The energy flux of a transmitted pulse at the distance \z\ = a/2, (a), (c) The non-resonant transmissiou by the slit 
(2a = 25 nm, b = 200 nm; t = 100 fs and 5 fs, respectively), (b), (d) The resonant transmission by the slit (2a = 25 nm, 
h = 350 nm; r = 100 fs and 5 fs, respectively). The central wavelength of the incident wave-packet is Ao = 800 nm. 




FIG. 5: The energy flux of a transmitted pulse at the distance \z\ — a. (a), (c) The non-resonant transmission by the slit 
(2a — 25 nm, b — 200 nm; r = 100 fs and 5 fs, respectively), (b) The resonant transmission by the slit (2a = 25 nm, 6 = 350 nm; 
T = 100 fs and 5 fs, respectively). The central wavelength of the incident wave-packet is Ao = 800 nm. 



